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INTRODUCTION 


Random  space  filling  problems  have  been  discussed  in  the  literature 
in  recent  years,  (for  example  [l],  [3],  [4],  [5],  and  [7]).  This  is 
an  interesting  and  usually  difficult  problem  in  geometrical  probability 
which  in  some  cases  is  motivated  by  specific  problems  in  physical 
chemistry,  biology,  and  other  disciplines.  Usually  one  is  concerned 
with  completely  covering  a fixed  region  by  random  regions.  However, 
the  specific  method  by  which  the  random  regions  are  placed  is  of  central 
importance  in  the  resolution  of  such  problems.  Because  of  the  complexity 
in  considering  this  in  any  general  fashion,  an  investigator  usually 
chooses  special  cases  to  examine.  In  this  paper  we  consider  some 
situations  considered  previously  by  other  authors,  expand  on  their 
analyses,  and  initiate  some  new  areas  for  which  some  solution  is  found. 

In  section  1 of  this  paper,  it  is  shown  that  the  moments  of  the 
random  variable  "proportion  of  the  fixed  region  covered  by  random 
regions"  form  a monotone  decreasing  sequence  that  converges  to  the 
probability  of  coverage.  This  is  shown  to  hold  in  great  generality 
and  exhibits  a link  between  the  method  of  Robbins  [5 ] for  calculating 
these  moments  and  the  random  space  filling  problem. 

In  Section  2,  formulas  are  developed  for  the  moments  of  the  propor- 
tion of  a set  that  is  covered  by  N independent  and  identically  distributed 
sets.  In  Section  3,  this  is  specialized  to  a particular  case  of 
coverage  of  a circle  by  random  arcs  of  random  size. 

Numerical  results  are  reported  in  Section  ^ for  probabilities  of 
coverage  of  the  circle  by  N uniformly  placed  random  arcs  of  random  sizes. 
Five  arc  length  distributions  are  considered  and  treated  by  Monte  Carlo 
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and  exact  methods. 


The  ordering  of  coverage  probability  by  decreasing  peakedness  of  arc 
length  distribution,  holding  expectation  constant,  is  explored  in  Section 
5.  This  conjectured  ordering  is  suggested  by  the  results  of  Section  4 
and  is  shown  to  hold  for  a large  class  of  arc  length  distributions  for 
which  exact  calculations  are  feasible. 
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1.  CONVERGENCE  OF  THE  MOMENTS  E Ak  TO  PROBABILITY  OF  COVERAGE 

Theorem  1 is  presented  here  in  full  mathematical  generality. 

For  ordinary  applications,  the  reader  may  wish  to  think  of  ^ as  a 
Euclidean  space  and  y as  Lebesque  measure  on  . It  is  in  this 
context  that  we  apply  the  theorem  in  subsequent  sections.  However^ 
the  full  generality  is  presented  here  because  the  theorem  was 
developed  that  way.  The  proof  of  theorem  1 is  essentially  the  same 
in  each  context . 

Let  u)  be  a measure  space.  That  is,  is  a a-algebra  of 

subsets  of  % and  y is  a countably  additive  set  function  mapping  ® to 
[0,°°].  Let  X be  a random  measurable  subset  of  such  that  its  measure 
y(x)  is  a random  variable.  Let  K be  a fixed  measurable  subset  of  * 
with  0 < y(K)  < 00  . Define  A = y(X0K)/y(K)  , the  proportion  of  K 
that  is  covered  by  X . Note  that  this  framework  includes  the  coverage 
problem  because  X may  be,  for  example,  the  union  of  N independent  and 
identically  distributed  regions.  To  avoid  measurability  problems,  we 
assume  y(XDK)  is  a random  variable.  The  event  "X  covers  K"  is  defined 
to  mean  y(XflK)  = y(K)  , that  is,  if  X covers  K up  to  a set  of 
measure  zero . The  main  result  linking  the  moments  of  A to  the  probabil- 
ity that  X covers  K is 

Theorem  1 : Each  moment  E A1”  (m  = 1,2,3, ...  ) is  an  upper  bound  on  the 

probability  that  X covers  K . The  sequence  of  moments  is  monotone 
decreasing, 

E A > E A2  > E A3  > •••  , 
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and  converges  to  the  probability  that  X covers  K : 


lim  E Am  P(X  covers  K)  . 

Moreover,  this  sequence  is  strictly  decreasing  except  in  tide  trivial 
case  ?(0<A<1)  = 0 - Thus  the  probability  of  coverage  is  determined 
in  a very  simple  way  from  the  moments  of  the  Proportion  of  K covered. 


Proof  of  Theorem  1:  Observe  that  A is  a random  variable  and  thus  a 

measurable  function  from  some  probability  space  (£2,^2, P)  to  the  interval 

[0,l]  and  that  P(A=l)  = P(X  covers  K).  0 < A < 1 implies  Am  >_ Am+1;  hence 
m m+1 

E A > E A and  the  sequence  of  moments  is  monotone  decreasing.  Next 
observe  that  for  each,  u e.  Q , lim  Am(oj)  = I^A  where  I^A  ^ 

is  the  indicator  function  of  the  set  {to:  ACw)=l}  . Moreover, 

| A (co ) | £ 1 . Thus  by  the  Lehesgue  dominated  convergence  theorem. 


lim  E A™  = lim  Am(w)dP(.u)) 

m->a>  m->oo 


I^A  l)(ta)dP(w)  = P(A=l)  • 

h 


Finally,  if  P(0<ACl)  f 0 then  P(Am  > Am+1 ) >0  and  EA”>E  Am+1 
for  each  positive  integer  m . □ 


As  an  example  to  show  that  measure-theoretic  problems  can  arise,  an 

example  is  now  given  in  which  y(x  f\  k)  is  not  measurable,  and  hence  is  not 
a random  variable. 

Let  - R,  the  real  numbers,  let  K = [0,l],  and  let  y be 
Lebesgue  measure  on  Let  E be  a nonmeasurable  subset  of  [ 0 , 1 ] . 
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For  a construction  of  E see  Royden  [6 ] p.  63.  The  random  set  X is 
defined  in  terms  of  a random  point  0)  selected  from  [0,l]  with  uniform 
measure.  We  set 


X 


|[o,i] 


[1,2] 


if  ti)  £ E 
if  0)  i E 


It  is  clear  that  the  value  of  X will  always  be  a measurable  subset 
of  and  that  y(x)  is  identically  1 so  that  y(x)  is  a (degenerate) 
random  variable.  However,  ^(X^Nk)  is  not  measurable  because 


{ 


uCxOk) 


if  0)  e E 
if  U)  £ E 


is  the  indicator  function  of  the  nonmeasurable  set  E.  This  is  presented 
here  to  highlight  the  measurability  aspect  referred  to  earlier  and  naturally 
this  situation  is  not  expected  to  occur  in  applications. 
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2. 


FORMULAS  FOR  THE  MOMENTS  IN  THE  CASE  OF  N I.I.D.  SETS 


Using  the  method  of  Robbins  [5]  the  moments  of  A = y(X O K)/y(K) , 
where  K is  a fixed  set  and  X a random  set  in  Rn,  y is  Lebesgue 
measure,  and  0 < y(K)  < °°,  are  given  by 


(2.1) 


E Am  = Iy(iO)]-m 


,x  e X ) dx_.  . . . dx 
* m 1 m 


To  avoid  measure— theoretic  pathologies,  we  assume  all  random  quantities 
we  write  are  random  variables.  This  will  be  guaranteed,  for  example, 
if  we  work,  in  the  framework  of  coverage  spaces.  A coverage  space 

defined  by  Ailam  in  [l],  is  a triple  <#>$>  M)  formed  by  two  probability 
spaces  d = (%  S,  P)  and  <B  = (^,  T,  Q)  together  with  a measurable 
set  Me  % x in  their  product  probability  space.  A random 

measurable  subset  M of  % is  determined  by  a random  point 

y 

y £ via  the  set  M by  taking  = (x  (x,y)  e M} . di  is  called 

the  coverable  space,  and  43  is  called  the  space  of  experiments.  A concrete 
illustration  that  is  instructive  in  viewing  this  concept  of  a coverage  space 
is  given  in  the  next  section. 

The  integrand  P ( x^ , ...  ^x^e  X)  is  difficult  to  evaluate  for  general  X. 


.N 


However,  in  the  case  X = X^  where  X^,...,Xjj  are  independently 

and  identically  distributed  sets,  the  integrand  may  be  expressed  in  terms 

of  x x and  X,  alone  as  follows. 

L m 1 
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Theorem  2 : 


m . N 

(2.2)  P(x. x eX)  = 1 + H-l)  l [P  (x  ,...,x  £X  )] 

k=l  J1  Jk 

N 

vhere  X = U ^ Xi  , and  X^-.-jX^  are  independent  and  identically- 
distributed  random  sets. 

From  this  we  immediately  obtain  a useful  formula  for  the  moments. 


Corollary: 

m , 

(2.3)  E Am  = 1 + l (-1)  (“)  [y(K)]"k 

k=l 

Thu  . to  find  the  mth  moment  of  A , we  need  only  calculate 

for  k = l,2,...,m  . Moreover,  after  finding  E Am  in  this  fashion, 

^ | ~| 

we  can  calculate  E A by  doing  just  one  more  integral,  namely 


1 IT 

[ p (.x^  > • • • ^ ) 3 * 

.Ik* 


[P(x 


1^- 


+1 


1’ 


II 


dx. 


.dx 


ia+1 


Proof  of  Theorem  2:  Let  I denote  the  indicator  function  of  the  set 

Y and  let  Y denote  the  complement  of  Y . Then 

P(x^,...,x^  e X)  = E Ix(x1)  ...  X^(x^)  • Since  X = , we  have 

IY(x  ) = 1 - II  I (x. ) . Thus 

x j i=l  xc  j 


P(x  , . . . ,x  e X) 
i.  m 


1 - 


N 

n 

i=l 


l 
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Expanding  the  first  product,  we  get 


m 

i + l (-1)1 

k=l 


N _ 

n 


1£J1< * • ♦ < i=l 


1 c(xj  > 
X.  J1 
1 


1 Ctxl 

XT  JI 

X 


Taking  expectation  inside  a finite  sum  and  using  independence,  we  get 

m , N r 

i + l (-i)k  l n e i rU  ) ...  i (x  )1  . 

k=l  l<j.  <m  i=l  L XT  X.  Jkj 

Theorem  2 now  follows  from  the  fact  that  each  X.  has  the  same  distribu— 

l 

tion  as  X^  and  that 


E I (x.  ) ...  I (x.  ) = P(x. 
X^  J1  X^  Jk  J: 


□ 
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3. 


SPECIAL  CASE:  ARCS  OF  RANDOM  SIZES  OH  A CIRCLE 


The  circle  has  held  a fascination  for  a number  of  authors  in 

connection  with  either  covering  the  circumference  with  random  arcs,  or 

covering  a region  with  random  circles.  There  is  an  interesting 

connection  between  the  two  problems  which  is  explored  here. 

Let  K denote  the  unit  circle.  Let  N arcs  X,  be 

1 N 

placed  on  K , with.  X__  centered  at  (h  radians,  and  of  length  L^  . 

We  assume  (^,...,0^  , L , . . . ,Ljj  are  independent  random  variables, 
with  each  Ch  uniformly' distributed  on  I0,2tt)  and  each  L^  distrib- 
uted with  a density-  h on  l0,°°)  . Let  A denote  the  proportion  of 

K covered  by  X = / , 2.  . We  wish  to  find  the  moments  of  A , and 
J i=l  i 

hence  P (X  covers  K.)  . 

This  problem  is  motivated  by  the  problem  of  covering  a plane  region 
with  random  unit  discs.  Observe  that  if  the  region  is  covered,  then 
each  disc  whose  closure  is  interior  to  the  region  must  have  its  boundary 
covered  by  other  discs.  Let  K denote  a fixed  unit  circle.  Place  K 
unit  discs  with  centers  uniformly  distributed  over  the  disc  of  radius- 
2 concentric  with  K . Hote  each  disc  determines-  a random,  sore  of  over- 
lap on  K (see  Figure  3.1).  It  can  be  seen  that  these  arcs  have  the 
distribution  mentioned  in  the  first  paragraph,  with  h(i)  = sin  l IjQ  ^jU). 
The  probability  that  the  H discs  cover  K Is  P(A=l)  . 

We  may  place  this  problem  within  the  framework  of  coverage  spaces 
as  follows.  Let  y be  the  surface  of  the  unit  sphere  in  3 dimensions 
parametrized  by  latitude  and  longitude  so  that 

= f(e,4>):  0 e [0,ir],  <f>  e [0,2tt)>.  Let  M = {(x,0,<j>)  e K x^: 
x e [$,0+$!  mod  2ir}.  Then,  as  outlined  in  section  2,  a point 
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(0,<j>)  e determines  a subset  M^g  ^ K.  In  this  case,  M^g  ^ 
is  the  arc  starting  at  <j)  and  extending  counterclockwise  with 
length  0.  If  we  place  uniform  measure  on  a straightforward 

calculation  shows  that  the  resulting  distribution  of  arc  length  is 
ho.  Figure  3.2  may  help  one  visualize  the  situation. 

First  we  calculate  the  moments  for  arbitrary  density  h_(J,)  . We 
will  use  Formula  (2.3)  as 


m 


(3.1) 


where 


E A™  = 1 + l (-1)K  © (2tt)  b 


k=l 


(3.2) 


b. 


[P(x1, . . . ^ Xx ) 3 dx1...dxfc 


Note  that  the  integrand  above  is  invariant  under  permutations  of 
(x^ , . . . ,x^  ) and  that 


[P(xl5. . . ,xk^X1)JN  dx1...dxfc_1 


K^l 


is  independent  of  the  value  of  x^  by  uniformity  of  C^.  . Thus  we 
may  set  x^  = 0 and  integrate  over  the  ordered  variables 

0 £X1  i i\  i i2,r  to 

(3.3)  bk  = (2tt ) (k-1 ) ! j (P^ , . . . ^^,0  dx.^ 

0<x^»**<xk  1<2tv 


Next,  transform  to  = x^  , a2  = x^x^ , . . . ,0^  ^ = x^_1  - x^g  > 


“k = ^ - \-i  ■ 2v  - Si  ai  • 


Note  that 


•*Vi 


ii 


(3.U)  P(x1,...,xk_1,0  ^Xx)  = £ P(XX  is  between  0 and  ou) 


i=l 


again  because  has  uniformly  distributed  center.  The  Jacobian  of 

the  transformation  ^ ) -*■  (a^,...,ak  is  1 so  (3.3) 

becomes 


(3.5)  bfe  = ( 2tt ) (k-1 ) ! 


r k 


R 


- N 

)■  P(X,  between  0 and  a.)  da,.. 

Li=l  1 1 


where  R = 


(a,  , . . . ,ol  , ) : all  a.  > 0 and  a.  < 2n  and 

1 k-1  l — ^i=l  i — 


rk.— 1 

= 2m  - • Next  we  calculate  P(X^  between  0 and  a)  for 

0 < a < 2m  . 


(3.6)  P(XX  between  0 and  a)  = E P(X^  between  0 and  ajL^)  = E 


(°-V+ 

2m 


where  (?)  + = max(^,0)  . has  density  h , and  we  integrate  by 

parts  to  obtain 


r a 


(3-T)  P(X_.  between  0 and  a)  = p— 

1 2m 


(a-£)hOt)d£  = 


J0 


2v 


a 


EUldl 


fZ 

where  H(ft)  = Jq  h(t)dt  is  the  cumulative  distribution  function  of  L 
Using  (3.7)  in  (3.5) 


(3.8)  bk  = (2m) (k-1)! 


2tt  l 
R i=l 


a. 

1 


IN 


HU)d& 


dai---dak-i 


For  completeness,  we  put  (3.8)  into  (3.1)  to  obtain 


m 


(3.9)  E Am  = 1+  l (-l)kQ(k-l)!  (27T)~(^k“l) 
k=l  k 


— I 

|2tt 

L i=l 


^i 


HU)d£, 


N 

da. 


13 


sin  £ 


Now  we  specialize  to  the  case  h(£)  = — X — Ir^  -i(£)  so  that 

o 2 L0,7tJ 


H(£)  = 
o 


1 - cos  £ 
2 

1 


0 < £ < it 

£ > tr 


and 


fa 


H(£)d£ 

J0  0 


a - sin  a 
Crr 

a 1_ 
2tt  ~ 4 


a tr 
a > 7T 


These  different  forms  on  [0,tr]  and  (tt,2tt]  complicate  matters,  hut 
difficulties  may  he  overcome  hy  breaking  the  region  of  integration  R 
up  into  subregions  where  some  (at  most  one!  ) a_^  is  greater  than  tr  , 
and  where  all  a^  are  smaller  than  tr  . Using  this  and  symmetry 
considerations.,  (3.81  becomes 


(3.10)  bk 


= C2TT)(k-^l)!  * k 


“ii.i  k 


ar-tr+tr  J2Cai-slji  ai! 

R n{ct1>ir} 


N 

da,,  ...  da.  , 
1 k-1 


+ (2tt  ) (k-1 ) ! 


■ 

k 

- 

(a „ - sin  a, ) 
i i 

1=1 

'R  H (all  a. 

<'tt} 

N 

da,  . . . da.  , . 
1 k-1 


k 

Recall  £ a.  = 2tt  , so  (3.10)  simplifies  to 
i=l  1 


l4 


(3.11) 


b = (2u)k! 
ic 


[ 


1 + 

if-  4tt 


_1_ 

4tt 


k 

l 

i=2 


1 N 


sm  a. 

l 


da. 


.da, 


k-1 


'R  fl{a1>Tr} 


+ (2TT)(k-l)l 


R 0 (all 


a <^jr} 

J 


For  k = 1 , we  may  use  (3.2)  directly  to  obtain 


(3.12) 


2tt 


For  k = 2 , we  use  (3-11)  and  the  fact  that  a^  and  = 2-rr-a-j^ 
cannot  both  be  less  than  IT  to  obtain 

r2TT 


b2  = Utt 


1 IT 


i+  _ 
4 4tt 


47  ^in(2iT— a^ ) 


Change  variables  and  simplify  to  obtain 


(3.13) 


b2  = 4 it 


r r 

1 

u* 

- w K + sin  “ij 

N 


da. 


This  may  easily  be  evaluated  numerically  for  any  particular  N . 

The  first  two  moments  of  A are  thus  (using  (.3.1)  5 (3.12)  and 

(3.13) ) 

3 N 

(3.14)  E A = 1 - {p  , 


K f77  r 

(3.15)  E A2  = 1 - 2*(J)  + M ■ 

3 1 1W 

7 - 47  (a  + sin  a)  da 
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The  third,  moment  is  more  of  a challenge. 


From  (.3 . 11 ) 


(3.16) 


For  ease  of  numerical  integration,  ve  change  variables  and  combine 
integrals: 


(3.17)  b3 


From  (3.1) 
(3.18) 


2 3 

Some  values  of  b^  b^,  E A,  EA  , E A and  variance  (A)  are 

2 

tabulated  in  Table  3.1.  Recall  from  theorem  1 that  E A,  E A , and 
3 

E A are  upper  bounds  for  the  probabilitv  of  coverage. 
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Table  3.1 


n_ 

E A 

2 

E A 

3 

E A° 

Var  (A) 

1 

it. 71 

22.67 

111.3 

.250 

.074 

.026 

.0118 

2 

3.53 

13.25 

52.1 

.438 

.211 

.111 

.0193 

3 

2.65 

7.89 

25.4 

.578 

.356 

.232 

.0218 

4 

1-99 

4.78 

12.9 

.684 

.488 

.361 

.0211 

5 

1 . 49 

2.96 

6.8 

.763 

.600 

.486 

.0186 

6 

1.12 

1.86 

3.7 

.822 

.691 

.592 

.0154 

7 

.84 

1.19 

2.1 

.867 

.763 

.681 

.0123 

8 

.63 

.77 

1.2 

.900 

.819 

.754 

.0096 

9 

.47 

.51 

• 71 

• 925 

.863 

.811 

.0073 

10 

.35 

.34 

.43 

.944 

.896 

.855 

.0055 

15 

.084 

.052 

.045 

.987 

.975 

.964 

.0011 

20 

.020 

.0092 

.0060 

.9968 

.9939 

.991 

.00022 

25 

.0047 

.0017 

.00093 

•99925 

.9985 

.9979 

.000044- 

IT 


COVERAGE  PROBABILITIES  FOR  THE  CIRCLE:  NUMERICAL  RESULTS 


Consider  again  the  problem  of  Section  3.  Let  F denote  the 
cumulative  distribution  function  of  the  arc  lengths,  and  let  P(n,F) 
denote  the  probability  that  n random  arcs,  with  uniformly  distributed 
centers,  cover  the  unit  circle.  In  this  section,  numerical  values  for 
P(n,F)  are  calculated  for  five  distributions  on  [0,tt]  , which,  like 
the  sine  distribution  of  Section  3,  are  symmetric  about  tt/2  . 

The  distributions  studied  are  as  follows: 


(l)  F , concentrated  at  tt/2 
/ ^ \ -rr,  . . . ...  sin  £ . 


(2)  FQ  , with  density  — - — Ij-  q ^jCL)  . Recall  that  upper  bounds  on 


PCnjF^)  were  found  in  Section  3. 

(3)  F^  , the  uniform  distribution  on  [0,tt]  . 


V*>  - 


1 /TT— A, 

1 ~ 


0 < £ < tt/2 


tt/2  < £ < tt  . 


This  was  chosen  as  an  intermediary  between  F^  and  F^  . 

(5)  F,.  , placing  mass  1/2  at  0 and  at  tt  . 

F^  through  F^  are  graphed  in  Figure  4.1. 

Exact  values  of  P(n,F)  were  obtained  for  F^  and  F^  . Computer 
simulations  gave  estimates  for  F^  , F^  and  F^  . 

P(n,F1)  is  found  using  Stevens  * formula  j"*'T'j  which  states  that  the 
probability  of  covering  a circle  with  n randomly  placed  arcs  of  fixed 
size  a is 


(4.1) 


[2ir/a]  . . n-1 

l (-nJ  (">  a - 

j=0  J ^ 
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ition  functions  of  arc  length. 


where  [t]  denotes  the  greatest  integer  contained  in  t . Thus 


(4.2)  P(n,Fx)  = 1 - n(^)n  1 + (g)^)11"1  - (^)(^)n_1 


P(n,F^)  is  found  by  conditioning  on  K , the  random  number  of  arcs  of 
size  tt  actually  placed,  and  using  (4.1).  Since  K has  the  binomial 
distribution  t$  (n,^-)  we  have 

P(n,F  ) = E P(n,F  |K)  = f (?)  2“n(l  - k(|)k_1)  . 

5 5 k=l  k 2 

To  simplify  this,,  note  that 


l k(£)(|)k_1  = 2(|-)n  l k(J)(i-)k(|)n  k 


k=l 


k=0 


k 3 ' 


2n  /3\n 
3 '“2.  * 


using  the  formula  for  binomial  expectation.  Thus 


(4.3)  F(n,F5)  = 1 - 2“n  (|-)n  . 

Monte  Carlo  results  were  done  on  Stanford's  IBM  370/l68  computer. 

For  each,  distribution,  10,000  simulations  were  performed  in  which 
random  arcs  were  placed  sequentially  until  either  the  circle  was 
covered  or  50  arcs  were  placed.  The  usual  unbiased  estimate  for  pro- 
bability of  success  in  Bernouilli  trials  was  used.  Numerical  values  of 
P(n,Fi>  are  tabulated  in  Table  4.1  and  graphed  in  Figure  4.2.  Note  that 
the  coverage  probabilities  vary  considerably  among  our  distributions,  even 
holding  number  of  arcs  and  expected  arc  length  constant.  For  example, 
with  6 arcs,  probability  of  coverage  ranges  from  .03  to  .27;  and  with  10 
arcs  it  ranges  from  .34  to  .62. 
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Table 

4.1.  Coverage 

probabilities 

for  n random  arcs  on  the 

circle. 

* 

t 

t 

t 

n_ 

P(n,F1) 

P(n,F2) 

P(n,F3) 

P(n,F^) 

P(n,F  ) 

1 

0.00 

0.00 

0.00 

0.00 

0.00 

2 

0.00 

0.00 

0.00 

0.00 

0,00 

3 

0.00 

0.001 

0.004 

0.02 

0.03 

4 

0.00 

0.01 

0.02 

0.06 

0.09 

5 

0.004 

0.03 

0.05 

0.12 

0.18 

6 

0.03 

0.07 

0.11 

0.21 

0.27 

7 

0.07 

0.13 

0.18 

0.29 

0.37 

8 

0.15 

0.21 

0.26 

0.38 

0.46 

9 

0.24 

0.30 

0.35 

0.46 

0.55 

10 

0.34 

0.39 

0.43 

0.54 

0.62 

11 

0.43 

0.47 

0.52 

0.61 

O.69 

12 

0.53 

0.56 

0.59 

0.68 

0.75 

13 

0.6l 

0.63 

0.66 

0.73 

0.79 

14 

0.68 

0.70 

0.72 

0.78 

0.83 

15 

0.74 

0.75 

0.77 

0.82 

0.87 

16 

0.79 

0.80 

0.8l 

0.86 

O.89 

17 

0.83 

0.84 

0.85 

0.88 

0.91 

18 

0.87 

0.87 

0.88 

0.91 

0.93 

19 

0.89 

0.90 

0.90 

0.93 

0.95 

20 

0.92 

0.92 

0.92 

0.94 

0.96 

21 

0.93 

0.94 

0.94 

0.95 

0.97 

22 

0.95 

0.95 

0.95 

0.96 

0.97 

23 

0.96 

0.96 

0.96 

0.97 

0.98 

24 

0.97 

0.97 

0.97 

0.97 

0.98 

25 

0.97 

0.98 

0.98 

0.98 

0.99 

* exact 
t simulated 
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5.  COVERAGE  PROBABILITIES  FOR  THE  CIRCLE  AND  PEAKEDNESS  OF  ARC  LENGTH 
DISTRIBUTION 


The  values  in  Table  4.1  and  Figure  4.2  suggest  that 
P(n,F1)  < P(n,F2)  <_  •••  <_P(n,F,_)  for  each  n . From  Figure  4.1  we  can 
see  that  the  F^  are  in  order  of  decreasing  peakedness  at  their  common 
expectation  tt/2  . We  say  F is  more  peaked  about  a than  G is 
about  3 if  Y^F  and  Z^G  imply  P(  | Y-a  | >_  t ) < P(  | Z-3|  ■>  t ) for 
every  t > 0 . This  definition  of  peakedness  is  due  to  Birnbaum  [2]. 

This  suggests  a conjecture,  namely  that  among  cumulative  distribu- 
tion functions  on  [0,2tt]  having  fixed  expectation  e , the  less  peaked 
F is  about  e , the  greater  P(n,F)  will  be  for  each  n . That- is, 

one  achieves,  greater  probability  of  coverage  if  one  chooses  a distri- 
bution of  arc  length  that  is  spread  out  rather  than  more  concentrated 
near  its  expectation. 

So  far  the  conjecture  is  supported  by  numerical  evidence  for  five 
distributions,  each  with  expectation  tt/2  . Next  we  study  a class  of 
distributions  that  lends  further  support  to  the  conjecture  and  enables 
us  to  prove  that  P(n,F1)  < P(n,F^)  . 

Let  F^  a denote  the  distribution  function  of  arc  length  L , 

where 


probability 


so  that  e = E L . For  this  class  of  distributions,  P(n,F)  may  be 
calculated  exactly  as  the  following  lemma  shows. 
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Lemma: 


(5.1) 


[2ir/a]  . 

P(n,Fe>a)  = 1 - Cl  - f)n  + I f-tf- 

j -*■ 


Proof : Conditioning  on  K , tlie  number  of  arcs  of  size  a placed, 

noting  that  K 'u($(n,— ) , and  using  Stevens’  result,  we  have 

Q. 


P(n,F  ) = E P(n,F  Ik) 

6,3,  6,3, 


n i i C 2ir/a  ] . , . , , 

- i oektL  - fr*  i (-d3 t^d  . 

k=l  k a a J=0  J 2Tr 


Separating  off  the  j = 0 term  and  rearranging  summations,  we  obtain 


Observing  that  the  inner  sum  is 


[l  - ;p-]n  E(^)  where  X 'U 
Tr  J 


J 


and  using  the  fact  that  X ^ ®(n,p)  implies  E(^)  = p^(a) 

J 0 

the  proof. 


completes 


□ 


We  will  be  interested  in  subclasses  of  distributions  with  expecta- 
tion e that  are  totally  ordered  by  peakedness  at  e . The  following 
proposition  exhibits  some  of  these. 

Proposition:  Let  e e (0,2tt)  . Then  the  class  = {Fg  e<_a<_2e} 

is  totally  ordered  by  peakedness  at  e . That  is,  e <_  a < b <_  2e 
implies  F is  more  peaked  at  e than  is  F , . 

6 , 3.  6 , U 


2k 


p |o 


Proof:  Let  Y 'u  F and  Z 'u  F where  e < a < b < 2e  . A 

0 9 a 6 j D — — - 

picture  will  help  one  visualize  the  situation  . 


We  must  show  P(|Y-e|  _>  t)  <_P(|z-e|  _>t)  for  all  t >_  0 . Consider 
three  cases: 

i)  if  t.  <_  b - e , then  P(|Y-e|  _>  t)  _<  1 = P(  |Z-e|  _>  t), 

ii ) if  b - e < t <_  e , then  P(  | Y-e  ( > t ) = 1 - — < 1 _ f- < p(  I z-e I > t ) , 

iii)  if  t > e , then  PC|Y-e|  _>  t ) = 0 = P ( | Z-e  [ _>  t).  D 

Theorem  3:  For  each  n , PCn.F^  &)  is  monotone  increasing  for 

a e [Tr/2,Tr]  . Thus  the  conjecture  holds  in  ^ . 

tt/2 

Proof  of  Theorem  3:  P(n,F^2  &)  assumes  different  forms  depending  upon 

whether  a e [tt/2,27t/3]  or  a e [2tt/3,tt]  . The  proof  given  here  is  for 
a e [2ir/35ir]  and  is  analogous  (but  more  complicated)  in  the  other  case. 

We  may  assume  n > 3 , for  if  n < 3 then  P(n,F  /0  ) = 0 when- 

IT  / ) Q. 

ever  a < tt  . Using  the  lemma. 


<5'2)  P(n’F„/2,a>  = 1 - - 


tt  -.n 
2a 


TTn 

2a 


\n-l  . tt n ( n-1 ) ,t\ 

V + 

e'  a 
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for  a e [27r/3,Tr]  . If  we  set  t = ir/2a  , we  need  only  show  that 


f(t)  = 1 - (l~t)n  - nt(|)n-1  + t(2t  - 1) 

2n 


1 3 

is  monotone  decreasing  for  t e [— , ^-]  . To  see  this,  first  observe 
that 


f”(t)  = n(n-l)[(|-)n~2  - (l-t)n_2]  >_0 


in  this  interval.  Hence  f’(t)  is  monotone  increasing  there. 

max 


I will  now  show  that 


rl  3,  f’(t)  = f'(f)  < 0 

t e 1 2 » 4]  U ~ 


Observe 


that 


<"<!>  - mzr  {<2)”‘1  + n - 1 - 'ff'1} 


where  we  define 


g(x)  = (-)X  + x - (|)X  . 


Observe  g(2)  = 0 , g(3)  = - r-  . For  x > 3 , 


g’  (x)  = (|)Xlog  \ + 1 - (|)Xlog  |ll-  (|)3log  | = -.37  < 0 


This  implies  f'(t)  < 0 , forcing  f(t)  to  be  monotone  decreasing  for 
t £ [|,|]  . 

The  proof  for  a e [ 7f/ 2 9 2rr/ 3 3 uses 


t D / -ri  \ -r  / 1 tTnII  to  ^ 3 .ii—l  tt  /n  wir  n ■,  tJL  /■  xi  \ / _\2 

P(n>FTT/2.a)  ‘ 1 - Cl-S>  + TO  ‘a1'?-11  - tr  (3>(' a -3) 


2 a 


2 6 


and  starts  by  observing  that  the  corresponding  third  derivative  is 
always  non-positive.  There  is  no  trouble  at  the  boundary  a = 2tt/3 
because  P(n,F  ) is  continuous  there.  ° 

TT  / 5 cl 

a)  is  plotted  in  Figures  5.1  and  5.2,  clearly  showing  its 
monotone  property.  In  figure  5«2  we  see  an  interesting  slope  dis- 
continuity at  arc  size  IT.  This  is  due  to  the  fact  that  two  arcs  of 
size  less  than  tt  cannot  possibly  cover  the  circle,  whereas  two  arcs 
of  size  greater  than  IT  can  cover  the  circle  with  positive  probability. 

Corollary:  PCnjF^)  < P(n,Fj_)  whenever  n >_  3 . 

Proof:  Observe  that  = F. v/2  vf2  and  F5  = Fir/2  it  * Hence 

PCnjF^)  <_P(n,Fc)  follows  immediately  from  Theorem  3 . Strict  inequality 

“I  O 

when  n >_  3 follows  from  f"(t)  > 0 for  t £ ^-]  in  the  proof  of 

Theorem  3 . D 

Theorem  Fix  e e (0,it)  . Then  P(.n,F  ) is  monotone  increasing 

..--n  r G 5 3. 

for  a e [max (e, it ) ,2tt]  for  each  n. 

\ * 

Proof  of  Theorem  4:  a > tt  , so  by  the  lemma  , 


(5.3) 


F(n,F  e,a)  = l-U-f)  (1-#) 


e sn-l 


Setting  t = — , 
a 


we  need  only  show 


is  monotone  decreasing  for  t e [gjr,  min(p  1)].  But 

f * (t ) = n jU-t)n_1  - (l-'lf)11-1]  1 0 
in  that  interval,  completing  the  proof. 

IT 

Corollary:  Let  — < e < ir.  Then  the  conjecture  holds  in  the  class 

of  distributions 

{F  : ir  < a < 2e} 

e,a  — — 

Proof : This  is  an  immediate  consequence  of  theorem  4 and  the  propo- 


sition of  this  section. 
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